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We study the Little-Parks effect for mesoscopic loops with very nonuniform thickness. The results 
follow the trend of the phase diagram obtained for almost uniform thickness. In particular, the 
singly-connected state is stable on a line segment delimited by two critical points. Most of this 
study considers loops with piecewise constant thickness; in this case the Euler- Lagrange equation 
can be integrated analytically. Under appropriate conditions, the temperature range where the 
singly-connected state is stable is proportional to the square of the ratio between the maximal and 
the minimal thicknesses. Our results may serve as a guide for planning experiments. 
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I. INTRODUCTION 

We deal with a loop of superconducting material as 
shown in Fig. |l| If the entire loop is superconducting 
and if it encloses a non-integer number of magnetic flux 
quanta, then single-valuedness of the order parameter 
requires the presence of some supercurrent /. This su- 
percurrent involves an energy price and, due to it, the 
transition to superconductivity occurs at lower temper- 
atures when the enclosed flux is p.on-integer. This effect 
was observed by Little and ParksEJ and an explanation in 
terms of the Ginzburg-Landau theory was provided by 
TinkhamE 

We have recently predicted that, if the thickness of the 
loop is not exactly uniform, then there exist situations for 
which superconductivity is broken at a layer, so that the 
superconducting part is actually singly-connectedB and 
no supercurrent flows. When this happens, we say that 
the sample is in the "singly-connected state" (SC). This 
is an interesting possibility, since it would allow for a new 
dimensionality of the regions where the order parameter 
may vanish. In the case of vortices, the order parameter 
vanishes along lines and there are claimsa that it cannot 
vanish on surfaces. On the other hand, it has been sug- 
gested that even for uniform thickness the SC state wr 
appear as an intermediate station in hysteretic paths. 
A systematic analytic study for the stability domain of 
the SC state in families of loops with thicknesses that 
deviate slightly from uniformity was carried on in Ref. ||. 
Mathematical justification for some of the assumptions 
in our model is given in Ref. [?]. Another situation in 
which the order parameter seems to vanish on a layer 
was considered in Refs. || - §. 

Since Ref. || may be too mathematically oriented, we 
shall review here its central findings. We assume that the 
thickness of the loop (radial dimension) is much smaller 
than the perimeter. In this case, for the range of temper- 
atures at which transitions occur, the order parameter is 
uniform on cross sections of the loop. (A "cross section" 
is a surface which is everywhere normal to the supercur- 
rent.) In addition, the free energy is insensitive to the 
magnetic field created by / and to the linear shape of 



FIG. 1. Thin loop of superconducting material with 
perimeter 2ttR and nonuniform cross section D. The magnetic 
flux enclosed by the loop is <E> and the supercurrent around it 
is /. 



the loop. We may therefore regard the loop as a circular 
ring with radius R equal to the perimeter divided by 2tt. 
We define by the angle which is obtained by deforming 
the loop to a circle, i.e. 6 is the arc length divided by R; 
the "thickness" D{0) is defined as the area of the cross 
section at 0. The origin = is defined by requiring that 
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The controllable physical coordinates are the tempera- 
ture T and the magnetic flux <3>. They enter our equations 
through 
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and 
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where £ is the coherence length, T c is the critical tem- 
perature in the absence of magnetic field, Tr is the tem- 
perature at which £ = R, $0 is the quantum of magnetic 
flux and n is some integer. 

Let ip(0) be the order parameter and ipo the order pa- 
rameter that would be obtained in the absence of mag- 
netic field. We define y (0) = \ip(6)/ipo\. The contribution 
of superconductivity to the free energy is proportional to 



j (-Ay 2 + -Ay 4 + (y ) 2 )Dd9 + (27rk) 2 A~ 1 , 



(5) 



where the prime denotes differentiation with respect to 9 
and 
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is a nonlocal term. (Some symbols, such as A and A^ 
are defined here not exactly as in Ref. §.) The term fc/A 
(cf. (||)) is proportional to the supercurrent I. Since (||) 
increases monotonically with /c 2 , the minimum of the free 
energy will always occur for k in the range [— ~, -]. Note 
that only the absolute value of ip(0) enters the expression 
for the the free energy; arg[^((9)] has been already worked 
out. 

A necessary condition for a minimum of the free energy 
is the Euler-Lagrange (EL) equation. For the SC state 
y(0) =0 and the EL equation is 



(Dy ) + XD(y - y 3 ) = 
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When y{0) is positive everywhere, we say that the sample 
is in the "doubly-connected state" (DC). In this case a 
better description is given in terms of w = y 2 and the EL 
equation reads 
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The phase diagram in the temperature - magnetic field 
plane is shown in Fig. ||. N, DC and SC denote the three 
possible states and Ti, Tn and Tin are critical lines at 
which second order phase transitions occur. When Tni 
is approched from below, y(0) decreases until it finally 
vanishes. Since Tm is located along |fc| = |, and the 




FIG. 2. Phase diagram. N: y = 0; DC: y > everywhere; 
SC: the order parameter vanishes at some place in the loop 
and becomes singly connected. Pi and P2 are critical points. 
A decreases with temperature and k 2 increases with the devi- 
ation from an integer number of flux quanta. 

most stable state is never at |fc| > the SC stability 
domain is restricted to a line segment. The ends of this 
segment are the critical points Pi and P2. When the 
magnetic field is varied and the line \k\ = ^ is crossed for 
small A (between Pi and P2), the current I vanishes and 
changes sign continuously. However, for large A (beyond 
P2), the current / changes sign discontinuously. 

Let us denote the positions of Pi,2 by (A = Ai 5 2, |fe| = 
2). Ai is given by the smallest value A for which the 
linearized Eq. (]?]) (without y 3 ) has nontrivial solutions. 
To locate A2 we define 
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Then, along Tin, Ar is negative (resp. zero, positive) 
when the value of A is less (resp. equal, larger) than A2. 
If D{0) is close to uniform, then, to first order in the 
deviation from uniformity, 

A 1 = i(l+/3); A 2 = i(l-2/3). (11) 

The analysis of Ref. |^ has an obvious shortcoming for 
experimental purposes: one would like the SC state to 
exist in a temperature range which is not too small, i.e. 
Pi and P2 are not too close. This occurs when \/3\ is not 
too small. But most of the results reviewed above were 
obtained under the assumption that D{0) is almost uni- 
form, which implies \/3\ << 1. The purpose of this article 
is to study forms of D{0) which are strongly nonuniform, 
to check whether the scenario depicted in Fig. || is still 
valid and whether only the first harmonic of D{6) af- 
fects the phase diagram [since D{6) enters Eq. (|ll| ) only 
through /3]. 



In most of this article (Sees. II and III) we study the 
case in which D(6) is a piecewise constant function. In 
this case Eq. (|8|) can be solved analytically. Sec. Ill con- 
siders the experimentally interesting limit of a weak link. 
We shall see that in this case P2 can move to tempera- 
tures where superconductivity is already well developed, 
provided that the ratio between the length and the width 
of the weak link is carefully tuned. In Sec. IV we consider 
deviations from the piecewise constant shape and Sec. V 
provides suggestions for experiments which, besides the 
appropriate form for D(9), are quite standard. 




II. PIECEWISE CONSTANT THICKNESS 

A piecewise constant thickness will be described by 
D(0) of the form 
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with < d < 1 and < #* < tt constants, and D{6) 
symmetric about = and about = tt. We are free 
to set the maximum of D(6) as 1, since the EL equation 
is invariant under multiplication of D(6) by a constant; 
d is the ratio between the minimum and the maximum 
thicknesses. 

In this case we have found a solution of Eq. (||) with 
the appropriate symmetry. For < 0* we require the 
order parameter to have a minimum at = 0. It has the 
form 

2v 2 

w{6) =A 1 - -^-raicn 2 (M,rai). (13) 
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For > 6* we require a maximum at = tt and the 
solution is 

2v 2 

w(9) =A 2 - — im 2 (l - m 2 )sd 2 (^ 2 6>, m 2 ). (14) 
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In these Repressions cn and sd are Jacobian elliptic 
functions V{ and mi are constants, and 
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With these forms, the constant term in Eq. (H) becomes 
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FIG. 3. Position of the critical points as functions of the 
ratio d = D(0)/D(tt) in Here 6>* = tt/2. 

Eq. (|la) for i = 1,2 and continuity of w and Dw' at 
= 0*H These equations are solved by Newton iterations. 
The integration in (||) is performed numerically. 

When w(0) = 0, A diverges and Eqs. ( |l6| ) simplify 
to A(2 - Ai) = 2vf. Still for w(0) = 0, the numerator 
and the denominator of the integrand in ( |lQ| ) vanish for 
6 — > 0. In this region, we expanded the integrand (ana- 
lytically) into a power series. 

It is sometimes useful to regard w(0) as an independent 
variable (and k as dependent). For < w(0) << 1, the 
integration in (|^) cannot be performed numerically near 
6 = 0. In this region, we substitute the integrand by a 
Pade approximant in powers of # 2 , with the denominator 
linear in 2 . 

As Ti is approached, mi j2 — > 0, and we recover the sit- 
uation in subsection 4.3 of Ref. |^. As A is increased away 
from Tj, mi and m 2 increase, but mi increases faster 
than m 2 . Not far from Ti, mi = 1. If A is increased fur- 
ther, mi and v\ become complex numbers, with \m\ \ = 1 
and arg(^i) = — ^arg(rai). Both regimes (mi and v\ 
real or complex) can be unified by writing mi = e~ 4t 
and v\ = re 1 : Eqs. ( p~3| ) and (|l5| ) together with relation- 
ships 16.9.1 and 16.11.3 of Ref. [l(J show that w is an 
even function of t. It follows that w(0) depends on t only 
through t 2 ; more precisely, w{6) is an analytic function of 
t 2 which is real when t 2 is real, and nothing special hap- 
pens when t 2 changes sign. This means that w is still real 
when mi and v\ become complex, and there is no singu- 
larity at mi = 1. Our computer programs used different 
parametrizations for the cases < m x < 1 and |mi| = 1. 
Since we are mostly interested in the region where A is 
large, in some cases we did not perform calculations in 
the region where < mi < 1. 



A. Mathematical Details 



B. Dependence of Ai,2 on d 



Eqs. (|13|)-(|16[) reduce the integro-differential equation 
(||) to the problem of determining the four constants V{ 
and rrii. We are thus left with four algebraic equations: 



By a straightforward generalization of subsection 4.3 
in Ref. pi we obtain 



d = tan( v^(tt - 0*)) tan( v / A7<9*) . (17) 

Ai is obtained by solving this equation. A2 is obtained by 
solving numerically Ar = for the situation w(0) = 0. 
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FIG. 4. Contour plot of k 2 for d = 0.2 and 0* = tt/2 near 
P 2 . /c 2 = \ along the line w(0) = 0. For A < A 2 = 0.858, /c 2 
decreases with w(0); for A > A2, k 2 "goes over a hill". 

Fig. H shows Ai and A2 as functions of d for 0* = it/ 2. 
The stability domain of the SC state increases monoton- 
ically with nonuniformity, following the trend given by 
the asymptotic formulae ([iT|). Fig. || shows a contour 
plot of k 2 as function of w(0) and A for 0* = tt/2 and 
d = 0.2. This plot looks as if had been copied from that 
obtained in Ref. for small nonuniformity. For A < A2, 
k 2 increases monotonically as w(0) decreases; the value 
\k\ = \ is reached when the order parameter describes 
the SC state, w(0) = 0. On the other hand, for A > A2, 
the value \k\ = \ is reached while w(0) is still positive; 
after that, for the minimum of the free energy, k will 
jump to — k when the magnetic field is swept, so that 
lower values of w(0) can only be achieved by metastable 
states. 



C. Dependence of A12 on 6* 

Fig. [5] shows Ai and A2 as functions of 6* for d = 0.1. 
The curve for Ai has a maximum deviation from j for 
0* = 7r/2, and is symmetric about this point. This is 
qualitatively what one would expect from (|TT|) . However, 
the curve for A2 has a marked maximum for << tt. 
This suggests that large domains of stability for the SC 
state may be found for weak links, provided that the 
dimensions of the link are properly tuned. The limit 

9* « 1 will be studied in Sec. III. 
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FIG. 5. Position of the critical points as functions 
of the fraction of the loop where D is small. Here 
D(0)/D(tt) = d = 0.1. The inset shows these positions 
for the shape given by Eq. (p5|), with p = 2. Again, 
D(0)/D(tt) « d = 0.1. 

D. Location of the Minimum of the Order Parameter 

According to Eq. (|l|), which was derived in Sec. 6 of 
Ref. ^, the angle where y(0) has a minimum and may 
eventualy vanish depends only on the geometry, and not 
on the temperature or the magnetic field. We want to 
check whether this result remains true when D(0) is far 
from uniform. The case described by Eq. (|l^ ) is not ap- 
propriate for this purpose, since the minimum will always 
be at = by symmetry. 

As a representative example, we consider a ring formed 
by three pieces of equal length, with cross sections in the 
ratio 2 : 10 : 5, as shown in the outer part of Fig. |[ We 
have calculated y(0;X) along the line Tni (see Fig. ||), 
where y vanishes for some 0. 

The inner part of Fig. ^ shows a polar contour plot 
of y(0; A). The radial coordinate is A2 — A, whereas the 
angular coordinate is just the angle along the ring. For 
A = Ai, y(0) vanishes for all (outer circumference); for 
A > Ai, y(Q; A) = along the nearly radial line on the 
upper right. As A increases, the layer where y{6) = 
moves towards the middle of the thin piece of the ring. 
This could be expected intuitively since, the smaller the 
coherence length, the less the influence of the neighboring 
pieces. The arrow at the right shows the angle where the 
order parameter would vanish according to the asymp- 
totic result ([j]). We see that the actual result is much 
less influenced by the neighboring pieces; this tells us 
that regions where D(0) is very small have a stronger 
influence than what is predicted if only their deviation 
from the average thickness is taken into account. 

We have also evaluated the order parameter along the 
line iy As we go from Pi towards A = k 2 = 0, the min- 
imum of the order parameter moves towards the middle 
of the thin piece, but only by a very small angle. 




FIG. 6. Outside: an asymmetric loop. Inside: contour plot 
of the order parameter in this loop for k 2 = - . 



III. WEAK LINKS AND SIMILARITY 

We want to evaluate the asymptotic value of A2 for 
0* << 1, near the region where A2 is maximum. We 
shall also assume (and verify a posteriori that this is a self 
consistent assumption) that A2 >> 1 and 0* /d remains 
finite in this regime. Our initial goal is the evaluation of 
Ar and, for it, we require w(0). w{0) is obtained from 
(|), with w(0) = 1/A = 0. Due to the continuity of Dw' ', 
iu'(#* + ) << w'(#* _ ). This boundary condition, together 
with A >> 1, make it reasonable to assume w{0) « 1 in 
the interval 0* < < n. In the interval < < 0* we 
write x = VXO, w{0) = w(x) = w(\f\0), and (||) becomes 



^w 2 + 2w 2 (l-w) =0, 



(18) 



where the derivatives are with respect to x and the 
boundary conditions are 

w(0) = 0, w{0) = 1, (19) 

where = V\0*. The similarity solution of 

w(x), is independent of d and depends on A and 0* only 

through the combination 0. 

By algebraic manipulation, we rewrite ( |To|) as 
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Noting that f(0) = \dw n \0)x 2 , that f(0) = 1 for > 0\ 
and changing the variable of integration, this becomes 
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At P 2 , Ar = 0. Using (fl|) and 0* « tt we obtain 

A 2 = 777777 N - / 7 v ^ . (22) 

v 7nD"(0) I (9 Jo iD(x)x 2 / v 7 



Eq. (|22|) gives a universal expression for the scaled 
as a function of only, for any d which is sufficiently 
small. Note that enters ( |22| ) not only through the limits 
of integration, but also through ([l9|). 
The solution of fll8| ) is of the form 
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where sn is a Jacobian elliptic function. Given m is 
obtained by numerical solution of ( jl9| ) and the right hand 
side of (g2|) can be evaluated. After dx/Xz is known, 
we also obtain 0*/d = 0/{dyf^). In this way we have 
obtained the universal curve for P 2 in Fig. 0, where 
has been swept as a parameter. We see that A 2 is of the 
order of 1/d 2 , provided that 0*/dis not too small. If this 
ratio is below ~ 7r, the assumption that w almost reaches 
1 at — 0* is no longer justified. 
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FIG. 7. Position of P 2 for the shape (|lg) with d « 1. 
Thick line: asymptotic line, given by (p2h ; in this case dy/X?, 
depends on 0* /d only. Thin line: d — 0.1. 

We are interested in large values of A 2 . Numerical 
maximization of © gives dxf\^ = 0.179 at = 0.883. 
Denoting by #* pt the value of 0* for which A 2 (d) is maxi- 
mum, this yields #* pt = l.blitd for d « 1. On the other 
hand, ( pTT| ) gives #* pt = 7r/2 for d « 1. A simple interpo- 
lation formula which has been found to give good results 
throughout is 



y opt 
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(24) 



We shall close this section by considering the behavior 
of the loop for "typical experimental conditions" near 



the SC state. For given d, we would like A2 — Ai to be as 
large as possible; for this purpose, a reasonable choice is 
#* = 6* pt . A typical temperature would be in the middle 
of the stability domain of the SC state. We therefore 
define A typ (d) = (Ai(0* pt ) + A 2 (<9* pt ))/2. Let us now fix 
#* = #* pt and A = A typ and sweep the magnetic flux 
(which determines k). As a response, the supercurrent 
I = k/A varies. It vanishes for k — and for k — |, and 
reaches its maximum / = J max for some intermediate 
value k = /c max . In order to characterize the behavior of 
the loop, we have evaluated k max (d). 

The criterion to determine k max (d) is as follows. Con- 
sider the system of equations that determines Vi and mi 
in Subsec. II. A. Note that k enters these equations only 
through the constant term in jl^ ) and that this term 
is an increasing function of /. Therefore, this constant 
term has a maximum possible value, which is attained at 
fc max . Above this value the system of equations has no 
solutions and, below it, it has two (one for k < /c max and 
one for k > k max ). At k = k max the system of equations 
has a double root and, therefore, its Jacobian vanishes. 
This provides the additional equation from which k max (d) 
can be determined. The results are shown in Fig. ||. Also 
shown are the values of k at which / has decreased by the 
factors y/2 and 2, as the SC state is approached. These 
are obtained by fixing the constant term in (16). 

Fig. H shows that the supercurrent decreases sharply 
as the SC state is approached. This effect might be used 
for calibration. 
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FIG. 8. Rate at which the supercurrent decreases when the 
SC state is approached under "typical experimental condi- 
tions" . Calculations have been performed in the region where 
mi in ( |l3| ) is complex. For d > 0.95 we have used the asymp- 
totic results developed in Ref. H. 



IV. DEVIATIONS FROM PIECEWISE 
CONSTANT THICKNESS 

Sections II and III answer our questions for the special 
case that D(6) is of the form (12). We ask now whether 
the behavior that has been found is typical of loops that 
resemble a weak link, or arises from spurious features of 



the model, such as the discontinuity at 6 = 6*. For this 
purpose, we consider D of the form 

D(6) = d + (1 - d) tanh p (<9/6>**) (25) 

for < 6 < 7r and symmetric about and tt. For #** << 
7r this shape models a weak link with D(0)/D(tt) = d. 
The larger the value of p, the more the shape of D will 
resemble a piecewise constant form. 



A. Smooth Thickness at = 

We have investigated the case p = 2 in Eq. (|2^). This 
time the differential equations in or (||) have to be 
solved numerically. Ai is found by linearizing ([?]), assign- 
ing an arbitrary normalization to 2/(0), and looking for 
the smallest value of A that gives a periodic solution. A2 
is found by requiring Ar = 0. 

In order to evaluate Ar, we have to compute y{6). For 
this we solve (0) with the boundary conditions y(0) = 
f/(7r) = 0. This is still unsufficient, since we need the 
ratio of the strongly vanishing terms f — f and // near 
= 0. For this purpose, a power series for y{6) was 
developed near 6 = 0. 

The inset in Fig. || shows Ai 5 2 as functions of #** for 
d = 0.1. Fig. ^| is an empirically scaled set of graphs 
for A2 (#**). The similarity obtained in Sec. Ill is no 
longer obeyed, but the qualitative results of the piecewise 
constant case are still true: for fixed d, A2(#**) has a 
maximum for some shape of the thickness, #** = #** t ; 
as d decreases, the peak of the unsealed A2(#**) becomes 
sharper, with larger A 2 (#** t ) and smaller #** t ; for #** < 
#** t , A 2 drops sharply. 
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FIG. 9. Position of P 2 for the shape (||), with p = 2. The 
lines are marked by the values of d. The scaling has been 
chosen empirically to mimic a nearly universal behavior, as in 
Fig. 7. 



B. D'(0) / 

We have investigated the case p = 1 in Eq. (|25|), which 
deviates stronger from piecewise constant than the case 
p = 2. The numerical analysis is very similar to that 
of the previous section, with one additional subtlety: 
strictly, the quantity Ar defined by (|^) and (|Io| ) should 
not be evaluated with the function w(6) which describes 
the SC state. Rather, we should evaluate it with w(0) > 
and take the limit w(0) — > 0. Let us define by Ar ? sc 
the value of Ar which is obtained by using the singly- 
connected order parameter in Eq. ([)]). Separating Ar 
into a part which is continuous when going to w(0) — > 
and a remainder, we obtain 



lim Ar 
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If D is smooth at = 0, the last term vanishes. 

U sing this formula, we have calculated the curves in 
Fig. |l0|. As might already be expected, these curves have 
maxima for suitable #**, and A2 increases when d de- 
creases, but the peaks are less pronounced and the de- 
pendence on d is weaker than in the previous examples. 
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FIG. 10. Like Fig. 9, for p = 1. 



V. DISCUSSION 

The examples considered here extend the results of 
Ref. |^ to cases where the thickness is far from uniform. 
The critical points Pi and P2 still exist and behave as 
in Ref. |[ The stability domain of the singly-connected 
state increases monotonically with nonuniformity. The 
dependence of Ai on 0* is symmetric about tt/2 and in 
qualitative agreement with Eq. (|ll|). 

There is a surprise in the dependence of A2 on 0*. For 
the thickness profile defined in (|l2|), \2(0) has a sharp 
peak when the ring has a "weak link" (<i, 6* << 1) and 
the ratio 9* /d is appropriately tuned, as found in Sec. III. 



In this case, the length of the temperature domain where 
the SC state is stable is inversely proportional to d 2 . 

Modern versions of the Little-Parks experiment are 
performed on mesoscopic loops, produced by microlito- 
graphic techniques. Smaller values of the perimeter 2ttR 
allow for smaller values of the coherence length and thus 
expand the relevant temperature range below T c . The 
basic mesoscopic loop is usually repeated a large number 
of times (e.g. Ref. [Ll]), but there are also cases in which a 
single loop is used (e.g. Ref. |l2|). The main contribution 
of this paper is the design of the thickness profile of the 
basic mesoscopic ring. 

Our goal is to design a sample with a large temper- 
ature domain for the SC state. Among the situations 
considered, we have seen that this is best achieved with 
the piecewise constant profile ([l2|), with #* = 0* pt . If #* pt 
is not accurately known, it is safer to take 6* > 6* pt than 
^opt- The temperature domain can be increased by 
decreasing the ratio d = D(0)/ D(tt). 

The piecewise constant thickness poses a mathemati- 
cal problem. Our one dimensional free energy in Eq. (||) 
is based on the assumption that the order parameter ip 
and the supercurrent density vector depend only on 
and remain constant on every given cross section. If 
D{6) is smooth, it has been proven that this assump- 
tion is approached when the thickness (linear dimension 
of D) is much smaller than the perimeter of the loop. 
But a piecewise constant thickness is not smooth; near 
the region where the thickness changes, the streamlines 
are strongly curved and the assumption that the current 
density remains constant on the cross section is an un- 
realistic picture. Nevertheless, the results of Sec. IV add 
confidence to the results of the piecewise constant case 
as a plausible limit. In order to have a safer answer to 
the piecewise constant case, and also since the usual ex- 
perimental situation is that the thickness is smaller than 
the perimeter only by a moderate ratio, a two dimen- 
sional treatment is desirable. Even if the results for the 
piecewise constant case turned out to stem from an over 
simplified free energy functional, the results of Sec. IV in- 
dicate that the temperature domain where the SC state 
is stable increases at least as the ratio D(tt)/ D(0). 

It is experimentally possible to make superconductig 
samples with weak links, such that D(0) is smaller 
than D(tt) by several orders of magnitude. However, in 
this case our theoretical treatment becomes inappropri- 
ate. The most important reason is that a large ratio 
D(tt)/D(0) implies a large volume where the streamlines 
are strongly curved and a large influence of the limita- 
tion discussed in the previous paragraph. Second, if the 
coherence length is much smaller than the perimeter, it 
may not necessarily be much larger than the thickness. 
And third, far from T c the Ginzburg-Landau theory is 
not necessarily a good approximation. Moreover, it is 
not necessary to have d « 1; a moderately weak link 
such as d ~ 0.1 will already give a temperature domain 
which is about twenty times larger than that of the "clas- 
sic" Little-Parks effect (0 < A < |), and should be easily 



accessible. 

Expression (||) for the free energy neglects the contri- 
bution of the magnetic field due to the supercurrent. This 
contribution was evaluated in Ref. ^ and it enhances the 
temperature domain for the SC state. This contribution 
may be important for a large periodic array of rings at 
A > A2 and \k\ = \\ in this case the ground state due to 
(||) is degenerate and a small perturbation may give rise 
to collective effects. 

Several possibilities for identification of the SC 
state andrnthe critical point P2 have been discussed 
previouslyM One possibility is to observe directly 
whether the order parameter vanishes at some layer 
6 = 0, by means of STM or some decoration technique. 
Another possibility is to check whether the supercurrent 
vanishes or jumps for \k\ « \\ this can be done by means 
of magnetometer, as in Ref. [l3|. The SC state was not 
observed in Ref. [H^; conceivably, the effective value of 
0*/d was too low. 

For A > A2 and \k\ = |, the free energy is a bistable po- 
tential. The potential barrier decreases and the overlap 
between the order parameters at the minima increases 
as A approaches A2. Therefore, fluctuations are expected 
between both minima, which occur more and more fre- 
quently as A — > Aj. However, for A < A2 the free energy 
as only one minimum, and no fluctuations are expected. 
This effect was indeed found in Ref. [H], but it was un- 
derstood that the sample is in the normal state after the 
fluctuations disappear. According to the present inter- 
pretation, at this stage the sample is in the SC state and 
is still superconducting. We have seen that P2 can be 
pushed to a regime where superconductivity is well de- 
veloped, so that this bistable potential might open the 
possibility for low dissipation logical devices. 

If P2 is pushed to a sufficiently low temperature and 
thermal excitations are rendered unimportant, one-may 
expect to achieve macroscopic quantum coherencelij of 
both minima in the bistable potential. This situation was 
looked for in the past with Josephson junctions .113 The 
setup described in this article seems to be cleaner and 
more promising: there are no uncontrollable parts in the 
loop, no external current and less sources of dissipation. 

To summarize: our results indicate that the tempera- 
ture region where the singly-connected state exists is not 
necessarily small, and can be readily found with existent 
experimental techniques. 
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